Balancing the Cosmic Energy Budget: The Cosmic X-ray 
Background, Blazars, and the Compton Thick AGN Fraction 



A. R. Draper and D. R. Ballantyne 

Center for Relativistic Astrophysics, School of Physics, Georgia Institute of Technology, 

Atlanta, GA 30332 

aden . draper@physics . gatech . edu 
ABSTRACT 

At energies > 2 keV, active galactic nuclei (AGN) are the source of the cos- 
mic X-ray background (CXB). For AGN population synthesis models to replicate 
the peak region of the CXB (~30 keV), a highly obscured and therefore nearly 
invisible class of AGN, known as Compton thick (CT) AGN, must be assumed 
to contribute nearly a third of the CXB. In order to constrain the CT fraction of 
AGN and the CT number density we consider several hard X-ray AGN luminos- 
ity functions and the contribution of blazars to the CXB. Following the unified 
scheme, the radio AGN luminosity function is relativistically beamed to create 
a radio blazar luminosity function. An average blazar spectral energy density 
model is created to transform radio luminosity to X-ray luminosity. We find 
the blazar contribution to the CXB to be 12% in the 0.5-2 keV band, 7.4% in 
the 2-10 keV band, 8.9% in the 15-55 keV band, and 100% in the MeV region. 
When blazars are included in CXB synthesis models, CT AGN are predicted to 
be roughly one-third of obscured AGN, in contrast to the prediction of one half 
if blazars are not considered. Our model implies a BL Lac X-ray duty cycle of 
~13%, consistent with the concept of intermittent jet activity in low power radio 
galaxies. 

Subject headings: galaxies: active — galaxies: jets — galaxies: quasars: general 
— X-rays: diffuse background 



1. Introduction 



Nearly half a century after the cosmic X-ray background (CXB) was discovered (IGiacconi et al 

19621 ). the majority of the CXB up to 10 keV has been resolved into distinct point sources by 
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20051 ). These discrete 



sources are active gala ctic nuclei (AGN), compact ex tra-galactic sources powered by accre- 
tion onto black holes (jLyndon-Belll Il969l ; iReea Il984j ) . As such, the CXB encapsulates the 
history of accretion onto super massive bla ck holes and provides a p owerful tool to aid scien- 
tific understanding of accretion processes (IFabian fc Barconslll992l ). It has been shown that 
a large portio n of this accretion is shrouded from our v i ew by intervening matter along the 



line of sight (ISetti fc Woltier 



1995: Fiore et al 



1999 



1989 



Celotti et al. 1992: Madau et al. 1994: Comastri et al. 



Treister et al.ll2009i ). For AGN spectral and spatial density mod- 
els to match the peak of the CXB at ~30 keV, the models must predict a larg e number of 
highly obscured sources known as Compton thick (C 
20031 : iTreister fc Urrvl I2OO5I : iBallantvne et alJbood 
hydrogen column density A^h ^1-5 x 10 24 cm" 2 (Treister et al.ll2009l ). making them practi 



) AGN (IRisaliti et al.lll999l : IUeda et al. 



Gilli et al 



cally invisible in the 2-10 keV band (jGhisellini et al.lll994j ). CXB synt hesis models predict 



2007). which have a neutral 



CT sources make up roughly half of t he obscured AGN population (IRisaliti et al.l Il999 



Ueda et al J 12003 : ITreister fc Urrvl 120051 : IBallantvne et al.ll2006l : iGilli et al.l 120071 ). 



In recent years, studies to observationally constrain the CT fraction have been under- 
taken. At first, small lo cal studies seemed t o agree with the model predictions that half of all 
ob scured AGN are CT (IRisaliti et al.lll999l : iGuainazzi et al.ll2005l ). However, a recent study 
by ITreister et al.l (120091 ). using samples from INTEGRAL and Swift observations and high 
redshift, IR-selected CT AGN candidates, suggests that CT AGN co ntribute only abou t 9% 
of the X-ray background, which contrasts sharply to the predicti o n by Gilli et al.l (120071 ) that 
CT AGN account for nearly a third of the CXB. iMalizia et all faoojl) studied 88 AGN ob- 
served by INTEGRAL/IBIS in the 20-40 keV band and found that at least 16% of obscured 
AGN are CT and >24% of the AGN in their local sample (z < 0.015) are CT. If ~75% of 
local AGN ar e obsc ured (IRisaliti et al.l 11999c ITreister et al.ll2009l ). the local AGN sample of 
Malizia et al.l (120091 ) suggests that >32% of obscured AGN are CT. Given the uncertainity 
of the CT AGN fraction, smaller classes of CXB contributors must be considered. The CXB 
contribution from the small class of AGN known as blazars has previously been ignored by 
CXB synthesis models and CT AGN fraction predictions, even though blazars are known to 
emit in a broad range from radio to TeV energies. To further constrain model predictions of 
the CT AGN fraction, the blazar class of AGN must be considered. 

Blazars are a unique and extreme class of AGN. Unified models of AGN, as summarized 
by Antonucci (1993) and Urry & Padovani (1995), explain blazars as radio galaxies with 
relativistic jets viewed close to the line of sight. Flat spectrum radio quasars (FSRQs) 
are relativistically beamed FRIIs (luminous radio galaxies) and BL Lac objects (BL Lacs) 
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are relativistically beamed FRIs (less luminous radio galaxies). The features which define 
blazars (extreme variability, high luminosity, high polarization, and radio core-dominance) 
are due to the relativi stic beaming caused by looking down the relativistic jet of the blazar 
( IPadovani et al.l 120071 ). The details of the spectral energy distrib ution (SED) of blazars is 
still a topic riddled with uncertainties (IKneiske fc Mannheimll2008l ). The extreme variability 
that distinguishes the blazar class ne cessitates simultaneous mult i- wavelength observations 
to understand the spectral properties (jGhisellini fc Tavecchioll2008l ). However, t he two-hum p 
form of the blazar SED, spanning from radio to 7-ray energies, is well known (jUrryl 119991 ) . 
The lower energy hump is due to synchrotron radiation w hile the higher energy hump is 



due to inverse Compton scattering (jUrry &: Padovani 



1995). It has been shown that blazars 



are significant progenitors of the 7 -ray background (jGiommi et al.ll2006l ; iNarumoto &: Totani 
20061 ; IKneiske fc Mannheiml 120081 ) ; therefore it is expected that blazars should have a non- 
negligible contribution to the CXB. 



Giommi et al.l (120061 ) predict that blazars should account for 11-12% of the soft CXB 
around 1 keV; however, no estimation is made for the blazar contr ibution to the peak region 
of the CXB around 30 keV. A recent study by lAjello et al.l (120091 ). based on the three year 
Swift/BAT blazar sample, claims that blazars contribute about 10 % of the X-ray background 
in the 2-10 keV band. In the 15-55 keV band lAjello et al.l (120091 ) predict blazars contribute 
~20% if blazars are modeled as a single p opulation or ~ 9 % if FSRQ s and BL L a cs are 
modeled as two distinct populations. Both iGiommi et al.l (120061 ) and lAjello et al.l (120091 ) 
found that blazars could contribute 100% of the CXB in the MeV band. 

Due to uncertainties in the low luminosity end of the AGN hard X-ray luminosity 
function (HXLF), multiple HXLFs must be considered (e.g., Ueda et al. 2003; La Franca 
et al. 2005; Silverman et al. 2008; Aird et al. 2009; Ebrero et al. 2009; Yencho et 
al. 2009) to understand the range of predicted CT AGN. Recent AGN HXLFs find that 
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minosity and density evolution (LADE), also fits the observational data well. Both the LDDE 
and LADE models are in keeping with the findings that the scarce, high-luminosity sources, 
quasars, show sharp posi tive evolution from z ~0-2, while less luminous sources, Seyferts , 
evolve more temperately (IBarger et al.ll2005l ; iBrandt fc Hasingerll2005l ; lHasinger et al.ll2005l ). 
Given the connection between AGN and galaxy evolution (e.g., Ferrarese & Merritt 2000; 
Smolcic 2009), it is not surprising that AGN evolution matches the trend of galaxy forma- 
tion, where massive galaxies formed earlier in cosmological time while smaller structures 
have waited until more recent times to form (e.g., Cowie et al. 1999). 
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In this work the blazar contribution to the CXB is predicted and the implications for 
the CT AGN fraction are discussed in the context of multiple HXLFs. In §2] we present 
the model used for the blazar and non-blazar AGN contributions to the CXB. In §0 our 
results are presented while discussions and conclusions are given in jjH We assume a A CDM 
cosmology with H = 70 km s" 1 Mpc" 1 , fi A = 0.7, and Vt m = 0.3 JSpergel et al.ll2007h . 



2. Calculations 



2.1. Blazar Contribution to X-ray Background 

2.1.1. Luminosity Function 

The widely accepted unified scheme of radio loud AGN is that FSRQs are FRIIs, lumi- 
nous radio galaxies, with jets pointed along the line of sight causing the observed radiation 
to be relativistically beamed , and that BL Lacs hay e an analogous relationship with FRIs, 
less luminous radio galaxies ( Urry Sz Padovani]|l995l ). Thus, the blazar luminosity function, 
which describes blazar space density and evolution, sh ould be well represented by the rel- 
ativistically beamed radio galaxy luminosity function (jUrry et al.l Il99ll ; iPadovani fc Urry 
19921 ; IPadovani et al.l[io07h . 



Willott et al.l (1200 ll ) represent the radio galaxy luminosity function as the sum of the 
low and high luminosity radio galaxies. The luminosity function is of the form 
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where L^imhz is the monochromatic luminosity at 151 MHz. 



Willott et al.l (120011 ) use three different evolutionary models for fh(z), which differ by 



the evolutionary scenerio for high luminosity radio galaxies for z > 2. Model A assumes a 
symmetric evolutionary scenerio where the density of high luminosity radio galaxies postively 
evolves until z ~ 2 and then evolves negatively at the same rate for z > 2. Mobel B assumes 
a postive evolution up to z ~ 2 and no evolution beyond z ~ 2. Model C assumes positive 
evolution of high luminosity radio galaxies up to z ~ 2 and for z > 2 negative evolution is 
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assumed, however the negative evolution is not assumed to be sym metric with respect t o 
z ~ 2. Due to the lack of high redshift sources in the sample used by lWillott et al.l (120011 ). 
all three models fit the data well. Model C is used here as it is the most general scenerio. 
The Willott luminosity function is converted from Einstein-de Sitter cosmology to ACDM 
cosmology. 

To relate the parent luminosity function to th e beamed lumino s ity fu nction we use 
the procedure laid out by lUrry fc Shaferl (119841 ) and lUrry fc Padovanil (I199ll ). The beamed 
luminosity L is assumed to be related to the rest frame luminosity C by 



L =(! + /<?)£, 



(4) 



where / is the fraction of the unbeamed luminosity that is relativistically beamed by the 
jet, and 5 = [7(1 — /? cos is the jet Doppler factor with (3 the apparent velocity in units 
of c, the speed of light. Finally, 7 = (1 — (3 2 ) 1 is the Lorentz factor, and 6 is the angle 
between the stream of the jet and the line of sight. The exponent p is defined as p = 3 + a for 
continuous jets and p = 2 + a for discrete jets, where a is the spectra l index, and account s 
for aberration, time contraction, and the blue-shifting of the photons (jlJrry Sz Shafer 1984 ). 
The probability of observing luminosity L given unbeamed luminosity C is 



P 7 (L|£) 



Pip \£ 



-(p+i)/p 



(5) 



The luminosity function for a particular Lorentz factor dQ^/ dlogL is found by integrating 
the intrinsic differential luminosity function d& e / dlogL such that 



d%{L) 
dlogL 



(6) 



To accommodate a range of Lorentz factors, we define the Lorentz factor distribution 71(7) 
so 71(7) oc 7 G and is normalized to one over the full range of Lorentz factors. Thus the 
observed blazar luminosity function dQ Q /d log L is 



d$ (L) 
dlogL 



d logL 



(7) 



By performing this procedure for L and £ pairings allowed by < 6 < 6 C , where 6 C is defined 
by f'5(6 c ) p = 1 (lUrry &: Padovanilll99ll ) the angle at which the beamed jet luminosity is equal 
to the intrinsic luminosity, we construct the luminosity function for the population of radio 
galaxies which have spectra dominated by relativistic beaming. 

To construct the blazar luminosity function the low luminosity function d<&i/d log L^imhz 
defined by equation [21 is set as the parent luminosity function for BL Lac objects and the 
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high luminosity function d<&h/d log Lisxmhz, denned by equation [3], is set as the parent lu- 
minosity function for FSRQs. Above Li 51MHz « 10 43 - 8 erg s _1 or Li 51MHz « 10 27 - 1 W Hz -1 , 
at z — 0, the luminosity function is dominated by FSRQs and be low this point th e BL Lac 
luminosity function is d ominant. For BL Lacs p = 3 + a = 2.7 (jUrry et al.lll99ll ). and for 
FSRQs p = 3 + a = 2.9 (IPadovani fc Urrylll992l ). To determine / and G, an a verage viewing 



angle and average Lorentz factor is selected in agreement with those found by iHovatta et al. 



(120091 ) . For BL Lacs the average viewing angle is 5.5° and average Lorentz factor is 10.3 and 
for FSRQs the average viewing angle is 4.4° and average Lor entz factor is 16.2. The range of 
Lorentz factors for both BL Lacs and FSRQs is 5 < 7 < 40 (IPadovani &: Urrvill992l) and the 
range of intrinsic luminosities (in erg s _1 ) is 39.5 < \ogLi 51MHz < 47.5 (jSmolcic et al.l 12009 



Ajello et al.ll2009l ). The results are not sensitive to the limits of integration as long as for 



the lower limit \ogL™™ MHz < 40.0 and for the upper limit ^ogL 1 [^ e MHz > 45.5. The FSRQ 
(solid line) and BL Lac (dashed line) luminosity functions at z — 1 are shown in Figure HJ 



2.1.2. Formalism 

We follow the formalism laid out in previous works (e.g., Comastri et al. 2005; Pompilio 
et al. 2000; Ballantyne et al. 2006) to compute the extragalactic background spectrum due 
to AGN and blazars and the associated number counts. The spectral intensity at X-ray or 
7-ray energy E due to blazars is given by 



c rw riogL TsTMHz d$(L 151MHz ,z) S E (L 151MHz ,z)d] 



2 



( )- ^jL Ll^ mh d\o gLl51MHz (i + znui + *Y + n k y/* dlogL ™ MH ^ 

(8) 

where d${Li§\MHzi z) / dlog Lisimhz is the blazar luminosity function at 151 MHz (in units 
of Mpc -3 ), Se(Li 5 imhz, z) is the observed spectrum at energy E (in units of keV cm -2 s _1 
keV -1 ) for a blazar with intrinsic 151 MHz luminosity Lx^xmhz at redshift z, and di is the 
luminosity distance of redshift z. A similar method is used to calcu late the number counts 



within a specified energy band as shown by iBallantyne et al.l (120061 ) equation 3 



2.1.3. Spectrum 



Blazars have a distinct spectral shape characterized by two bumps, a synchrotron peak 
and a higher frequency inverse Compton (IC) peak. The 'Blazar Sequence', in which the 
spectrum of a bla zar can be uniquely determined based solely on the bolometric lum inos- 



ity of the blazar (IFossati et al.l Il997l Il998l : iGhisellini et al.l Il998l : iDonato et al.l 12001( 1 was 



considered as a blazar SED model; however, the anti-correlation between the synchrotron 
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peak frequency a nd the synchrotron peak luminosity has been proposed to be due to sample 
selection effects ( jCaccianiga &: Marcha 2004; Pa dovanil 120071 ; iGhisellini fc Tavecchiol 120081 ) 
and possibly variable Doppler boosting (INieppola et al.ll2008l ). A 'new Blazar Sequence' has 
been prop osed which parametrizes the blazar SED by the black hole mass and the accre- 
tion rate (iGhisellini fc Tavecchiol 120081 ) ; however, the new Blazar Sequence still predicts an 
anti-correlation betwe en the synchrotron peak frequency and synchrotron peak luminosity 
(INieppola et~aJl2008h . 



As there is no widely accepted blazar SED model, the blazar population was split into 
three subclasses (FSRQs, Low-peaked BL Lacs (LBL), and High-peaked BL Lacs (HBL)) 
and SED models were created for each subclass based on well sampled spectra of subclass 
members. Details of the SED parametrization can be found in Appendix |A] For each 
subclass we set the synchrotron peak frequency, us, the IC peak frequency, uic, and the 
relative luminosity of the synchrotron and IC peaks, Lic/Ls- 

For FSRQs the sy nchrotron peak frequency is set based on the average peak frequency 
depicted in Figure 5 of lAnton fc Browne! (120051 ) for the 1 Jy sample. This gives an average 
FSRQ peak freque ncy of 10 14 Hz. The IC peak frequency of FSRQs is set at 1 MeV 



in accordance with lAjello et al.l (120091 ) . To determine the ratio Ljc/Ls the literature was 
searched yielding the spectra of forty FSRQs shown in Table HJ In cases where multiple 
spectra are given for the same object, the spectra corresponding to the most quiescent state 
is considered. Most sources show ratios of \og(L IC / L s ) ~ 1.0 to 2.0 with a few sources 
showing ratios as low as \og(Ljc/ Ls) ~ 0.0 or as high as log(Lic/Ls) ~ 3.0. The mode 
value for the ratio log(Lic/Ls) = 1.0 is used. 

The BL Lac obj ects are divided into HBLs and LBLs assuming 10% of the BL Lac 
population are HBLs (jPadovani et al.ll2007l ). Similarly to the FSRQ subclass, the HBL and 
LBL synchrotron peak frequencies were assigned the average peak frequency from the Slew 
sample and 200mJy sample, respectively, as depicted in Figure 5 of lAnton fc Brownd (120051 ). 
The average peak frequency is 10 16 5 Hz for HBLs and 10 14,5 Hz for LBLs. To set the IC 
peak frequency and IC peak luminosity for the HBL and LBL spectra,the ratios vicl v s 
and Ljc/Ls from the SEDs modeled in the literature were consulted (see Table [2]). All 
sources had \og{i>ic /vs) ~ 8.0 to 9.0 and \og(L IC / L s ) ~ -1.5 to 1.5. The mode values 
for the ratios were used, giving log(z// cl v s)= 8.0 and Xo gjLjcl Ls)=§£$. Tw o LBLs were 
found to have complete spectrum models, BL Lacertae (IBerger et al.l 120081 ) and 3C 66A 
( jJoshi fc Bottcherl 120071 ). The spectra of BL Lacertae and 3C 66A both have logiyic l^s) ~ 
7.0 and log{L IC /L s ) « 0.0. 



Figure [2] depicts the spectra used for FSRQs (solid lines) for Li 5 imhz 



10 43.o erg s -i ; 



LBLs (dashed lines) for Li 51 mhz= 10 erg s , and HBLs (dot-dashed lines) for L\§\mhz 
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10 40.0 erg s -l_ 

2.2. AGN Contribution to the CXB 

The non-blazar AGN contribution to the hard CXB is computed using standard syn- 
thesis modeling techniques (e.g., Comastri et al. 1995; Treister & Urry 2005; Ballantyne et 
al. 2006; Gilli et al. 2007). The fraction of Type 2 AGNs, f 2 , is assumed to be a function 
of both redshift and 2-10 keV luminosity, L x : / 2 oc (1 + z) a (\ogL x )~ b , with a = 0.4 (Bal- 
lantyne et al. 2006; Treister & Urry 2006) and b = 4.7. This evolution is normalized so that 
the non-blazar AGN type 2 to type 1 ratio is 4:1 at z = and logLx = 41.5. The redshift 
evolution is halted at z — 1 in analogy with the evolution of the cosmic star-formation rate 
density (e.g., Ghandi & Fabian 2003; Hopkins & Beacom 2006). 

AGNs with absorbing column densities less than log Ah = 22 are considered to be 
unabsorbed type 1 sources, and they are distributed evenly over the following columns: 
log Ah = 20,20.5,21 and 21.5. Compton-thin Type 2 AGNs are also distributed equally 
over log A^h = 22, 22.5, 23 and 23.5. As current hard X-ray luminosity functions are missing 
Compton-thick AGNs, a parameter controlling the Compton thick fraction is defined, and 
any CT AGNs are distributed equally over log A H = 24, 24.5 and 25. It is assumed that CT 
AGNs evolve in the same manner as less absorbed AGNs. The CT fraction, f CT is defined 
to be the fraction of obscured AGN which are CT. 

The unabsorbed rest-frame AGN spectrum consists of a power-law with photon index 
T and an exponential cutoff at energy E cut , combined with a neutral reflection component 
calculated using the 'reflion' model within XSPEC (Ross & Fabian 2005). The strength 
of the reflection features in the total spectrum is typically parameterized by a reflection 
fraction, R, that is related to the covering factor of the reflector. Observationally, it is found 
that the strength of the reflection features decreases with luminosity, a relationship that 
is sometimes called the X-ray Baldwin effect (e.g., Bianchi et al. 2007). Therefore, we do 
not assume a constant value of R, rather the reflection spectrum is added to the power-law 
component such that the equivalent width of the Fe Ka line agrees with the observed X-ray 
Baldwin effect found by Bianchi et al. (2007). In this way, the observed decrease in the 
strength of the reflection features with luminosity can be naturally included in the synthesis 
model. The relationship by Bianchi et al. (2007) gives an Fe Ka equivalent width of 143 eV 
at log Lx = 41.5, which is approximately equivalent to R = 1.1 for a T = 1.9 spectrum. The 
reflection fraction is proportional to the Fe Ka EW, so R reduces to approximately 0.4 at 
log Lx = 44 and 0.1 at log Lx = 47. All the models presented here assume E cnt = 250 keV. 
Finally, following Gilli et al. (2007), spectra with T = 1.5 up to 2.3 are Gaussian averaged 
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around T = 1.9 to account for the observed dispersion in AGN spectral slopes. This results 
in a final rest-frame spectrum with luminosity Lx with the correct reflection strength. 



2.3. AGN HXLFs 



The AGN HXLFs considered are listed with their parameters in Ta ble El Five of the 



HXLF s considered are LDDE models as proposed by lUeda et al.l (120031 ). The lAird et al. 



(120091 ) HXLF considered is a LADE model. For the LDDE models, the luminosity function 
is described by the local luminosity function, d$(Lx, z = 0) / dlog Lx and an evolution factor 
e(z, L x ), such that 

d$(L x ,z) d$(L x ,0) 



d log Lx d log Lx 
The local luminosity function is of the form 



-e(z, L 



Xj 
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d log L x 

And the evolution factor is given by 
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L x > L a 
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(9) 



(10) 



(12) 



The LADE model takes the same form as the LDDE local luminosity function defined in 
equation [10] and and A are allowed to evolve with redshift such that 

r/1.0 • :,Y' : . /1.0 • :, x 

logL^^;) = logLo - log 



and 



1.0 + z 



log A(z) = logA + a(l + z). 



1.0 + 2 



(13) 
(14) 



3. Results 

3.1. Blazar contribution to the CXB 

Figure [3] shows the calculated blazar and AGN contribution to the CXB and 7-ray 
background. The dotted line shows the AGN contribution to the CXB assuming the HXLF 
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presented by lUeda et aJj (120031 ) and far = 0.5. The dashed line is the FSRQ contribution 
while the dot-dashed line shows the BL Lac contribution. The BL Lacs are found to con- 
tribute greater than 100% of the 7-ray background. The possibi lity that the oyeresti mation 
was due to inappropriate beaming parameters was investigated. iHovatta et all (120091 ) found 
that BL Lacs tend to have Lorentz factors of 1.0 < 7 < 38. When this range of Lorentz 
factors is used and the BL Lac beaming parameters are set so that the average Lo r entz f actor 
is 10.3 and the average viewing angle is 5.3°, in agreement with IHovatta et al.l (120091 ). the 
model overpredicts the 7-ray background by a factor of ~200. 

We have assumed that all radio galaxies viewed in the appropriate orientation are 
blazars. However, evidence suggests that AGNs, specifically radio galaxies, are an inter- 



mittent phenomenon ( 


Burns et al. 


1983 




Cavaliere & Padovani 


1989 




iloettieer et al. 


1993; 


Franceschini et al. 


199a 


Schoenmakers et al. 


1999. 


200(1 


Venturi et al. 


2004; 


Jamrozv et al. 


2007; 


Parma et al. 


2007) 


. BL Lac high energy radiation is primarily produced through syn- 



chrotron self-Compton (SSC) upscattering, since BL Lac jets propagate through regions with 
very little external radiation (iGhisellini et al.ll2009al ). For BL Lacs to emit high energy ra- 
diation, processes within the jet must accelerate electrons to relativistic speeds such that 
the electrons have enough energy to create synchrotron photons and then upscatter those 
photons. Therefore the high energy IC component of BL Lac jets would perhaps only be 
significant du ring infrequent event s that cause rapid acceleration of electrons to high ener- 
gies. Indeed, iGiommi et al.l (120061 ) found that unless BL Lacs have a small duty cycle the 
predicted blazar 7-ray emission would over predict the 7-ray background. The blazar and 
AGN contribution to the CXB assuming an X-ray duty cycle of 13% for BL Lacs, as shown in 
Figure HI fits the data well. The BL Lac number counts predicted by this model in the 15-55 



keV b and are shown in Figure [5] with the BL Lac number counts observed by lAjello et al. 



( 120091 ). The number coun ts predicted b y this work are slightly larger but within a factor of 
2.5 of the observations by lAjello et al.l (120091 ) . 



More powerful sources, like FSRQs, are belie ve d to accrete more efficiently and at a 
higher rate than low luminosity sources (IHol 120081 ) . iBauer et al.l ( 120091 ) studied variability 
in Palomar-QUEST survey blazars and found evidence that FSRQ duty cycles are greater 
than BL Lac duty cycles. No duty cycle for FSRQs is accounted for h ere, and ind e ed th e 
background due to FSRQs is in reasonable agreement with that found by lAjello et al.l (120091 ). 



In the soft X-ray band (0.5-2 keV) t he blazar cont r ibutio n is found to be ~12%, in 
agreement with the prediction of 11-12% by lGiommi et al.l (120061 ) . The blazar contribution in 
the hard X- ray band (2-1 keV) is found to be ~7.4%, in rough agreement with the prediction 
of ~10% bv lAiello et aD (|2009[ ). In the 15.0-55.0 keV range blazars contribute ~8.9% of the 
X-ray background, in good agreement with the prediction of two distinct blazar classes by 
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Ajello et al.l (120091 ) of ~9%. Emissions from B L Lacs is found to account for the MeV 



background, in agreement with previous works (iGiommi et al.l 120061 ; iNarumoto &: Totani 



20061 : iKneiske fc Mannheim! 120081 ). 



3.2. Implications for CT AGN 



When the contribution of blazars to the CXB is properly consider ed, fewer CT AGN 
are required. In Figure H] the AGN contribution to the CXB, as given by lUeda et al.l (120031 ). 
is shown with far = 0.4, in contrast to the canonical fcr = 0.5. The CT fraction, fcT 
required to appropriately model the peak of the CXB is shown for various HXLFs in Table 
H]as well as the for needed if the contribution of blazars is not considered. There is a ~10% 
uncertainty in the peak intensity of the CXB (e.g. HEAO-1 vs Swift). For the purposes of 



this work we assume EFe ~ 44.2 keV cm 2 s 



- 1 str" 1 



at 30 keV. Note that the estimated 



CT AGN fraction strongly depends on the CXB peak intensity adopted in the model. The 
Yencho and Silverman HXLFs depend heavily on CT sources to match the peak of the CXB; 
therefore, when blazars are considered the CT AGN fraction is still quite high. The number 
density of CT AGN at z = with Lx > 10 43 erg s -1 is also shown in Table H] for the case 
of no blazar contribution to the CXB and for the case of blazar contribution to the CXB as 
described here. The number density as a function of redshift of CT AGN with Lx > 1 43 erg 



s neede d to model the peak of the CXB for the HXLF given by lUeda et al. (I2003h (solid 
lines) and lEbrero et al.l (120091 ) (dashed lines) are shown in Figure El The thin black lines are 
the case where the blazar contribution to the CXB is not taken into account. The thick blue 
lines show the case where blazars are considered. The CT AGN density is reduced by a factor 
of 1.7 for the Ueda HXLF and a factor of 3.0 for the Ebrero HXLF. The AGN luminosity 
functions proposed by Silverman et al. (2008) and Yencho et al. (2009) require CT fractions 
where the majority of obscured AGN are CT. The AGN luminosity functions proposed by 
Ueda et al. (2003), La Franca et al. (2005), and Ebrero et al. (2009 ) re quire less than 



half th e obscured AGN are CT, in agreement with iMalizia et al.l (120091 ) and iTreister et al. 
(l2009h . 



4. Discussion and Summary 

It is clear that blazars make a non-negligible contribution to the CXB and significantly 
reduce the number of CT AGN predicted, and may be primarily responsible for the MeV 
background. This paper presents an upper limit to blazar contribution to the CXB by 
utilizing the unified model of radio-loud AGN. The main conclusions found here do not 
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change with a different choice of AGN radio luminosity function (e.g., Condon et al. 2002; 
Sadler et al. 2002; Best et al. 2005; Kaiser & Best 2007; Mauch & Sadler 2007); however, 
beaming parameters need to be modified as these luminosity functions do not distinguish 
between FRIs and FRIIs or high and low luminosity sources. 



A recent study by lAjello et al.l (120091 ). using the three year sample of Swift/BAT blazars, 
finds similar results for FSRQs as those found here; however, this work finds a greater 
contribution to the CXB and cosmic 7-ray background by BL Lacs. Due to the small num ber 
statistics and small redshift range of the Swift/BAT BL Lac sample, lAjello et al.l (120091 ) are 
not able to uniquely determine the evolutionary parameters and thus assume no evolution 
for BL Lacs. This work assumes BL Lacs evolve in the same manner as low luminosity 
radio galaxies. Also, lAjello et al.l (120091 ) assume a simple power law SED model for BL Lacs 
whereas this work utilizes an SED model based on average BL Lac properties taking into 
account the variety of BL Lac subclasses of LBLs and HBLs. Severa l studies have found 
that BL Lacs contribute substantially to the cosmic 7- ray background (jGiommi et al.l 12006 



Narumoto fc Totanill2006l ; iKneiske fc Mannheiml 12008), thus it is not ex pected that the BL 



Lac contribution to the CXB is negligible, as found by lAjello et al.l (120091 ). As this work uses 
a more physical BL Lac SED model and a reasonable evolutionary model, we expect that 
this work may more accurately model the BL Lac contribution to the CXB, although the 
factor 2.5 discrepancy in the BL Lac source counts in the hard X-ray band must be solved 
in future works. 



Treister et all (I2009h find the density of CT AGN at z = wit h L x > 10 43 erg s" 1 
is ~2.2 x 10~ 6 Mpc~ 3 . The luminosity function of lUeda et al.l (120031 ) predicts the density 



of CT AGN with L x > 10 43 erg s" 1 at z = to be 7.3 x 10~ 6 Mpc~ 3 if blazars are not 
considered and 4.4 x 10 -6 Mpc~ 3 if bla z ars ar e considered. With the blazar contribution 
to the CXB cons i dered , the lUeda et al.l (120031 ) over predicts the CT AGN density found 

byT 



Treister et al. 



Ebrero et al 



_ (120091 ). by a factor of 2. Conversely, the luminosity function proposed by 
J2OO9J) predicts the density of CT AGN with L x > 10 43 erg s" 1 at z = to 
be 1.1 x 10~ 6 Mpc -3 if blazars are not considered and 3.6 x 10~ 7 Mpc~ 3 if bl azars are 



considered, which is a factor of 6 smalle r than the den s ity re ported by lTreister et al.l (120091 ) . 
According to the INTEGRAL results of iMalizia et~aD £53), the f CT > 0.32 with no upper 
limit given. Between different HXLFs there is a large scatter in the predicted far and 
the predicted CT AGN density varies by a factor of 30. This clearly illustrates the limits 
imposed by the uncertainty of the low luminosity end of the AGN HXLF and how important 
it is for future missions to probe this portion of the HXLF. 

It has been shown that blazars, specifically BL Lacs, contribute the majority of the 7 - 
ray background (jGiommi et al.ll2006l ; iNarumoto &: Totanill2006l ; IKneiske fc Mannheimll2008l ) . 
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Giommi et al.l (120061 ) found that unless BL Lacs have a small high energy duty cycle the 
predicted blazar 7-ray emis sion would over predict the 7-ray background. Furthermore, 
Kneiske fc Mannheiml (120081 ) suggests that using radio blazar luminosity functions may cause 
an overestimation of the number of sources emitting robustly at higher energies, as it is not 
certain that all radio sources will have strong X-ray and 7-ray emission. Physical and 
evolu tionary mod e ls of quasars indicate that AGN activity is short-lived and possibly recur 



rent (Soltanlll982l ; ICavaliere fc Padovanilll989l ; IChokshi fc Turner! 119921 ). iFranceschini et al. 



(119981 ) showed that long-lived, continuous AGN activity is not consistent with the black hole 



mass function they calculated from their sample of 13 local galaxies, but short-lived and 



AGNs have been observed (IBurns et a . 


1983: 


Roettieer et al. 


1993 




Schoenmakers et al. 


1999; 


Venturi et al. 


2004; 


Jamrozv et al. 


2007 




Fabian et ali 


2009). Sources have also been 



observ ed which have relic radio lobes but the AGN activity is not curr ently in an active 
phase (IParma et al.l 120071 ; iDwarakanath fc Kald |2009| ; iFabian et al.l 120091 ) . Recent observa- 
tions by Hubble Space Telescope and Chandra of the relativistic jet of nearby M87 provide 
evidence for the intermittent nature of jet X-ray emission (e.g., Perlman et al. 2003; Har- 
ris et al. 2006; Stawarz et al. 2006; Madrid 2009). Large amplitude flaring h as been 
obser ved from the previously quiescent knot HST-1 in the jet of M87 since 2000 (Madrid 
20091 ). This flaring activity is shown to be consistent with shocks occurring within the jet 
as faster mo ving particles collide with slower relativistic particles inj e cted into the je t at an 
earli er time JPerlman et ali l2003l: Istawarz et aDbood : iMadridl 120091 ) . IPerlman et all J2003f ) 
and IStawarz et al.l (120061 ) suggest the recent X-ray flaring of HST-1 is directly related to 
material injected at the base of the jet 30-40 years ago. Therefore, a jet X-ray duty cycle is 
expected. 

Finally, due to the spectral steepening that occurs after the flow of energetic particles 
into the jet has ceased, the b est frequency rang e to search for relic radio lobes is the low 



radio regime, less than 1 GH z (IParma et a 



luminosity function given by IWillott et al 



20071) . Therefore, it is likely that the radio AGN 



(120011 ) at 151 MHz includes relic radio lobes. As 



this would affect the low luminosity end of the luminos ity function more prevalentl y as relic 
sources tend to not be as luminous as active sources (IDwarakanath fc Kald 120091 ). the BL 
Lac luminosity function found here may overpredict the number of BL Lacs. Thus, the 
average BL Lac X-ray duty cycle is likely to be somewhat larger than the 13% found here. 



The authors thank C.M. Pierce for reviewing a draft of this paper and the referee for 
helpful comments. 
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A. Blazar Spectrum Parametrization 



The synchrotron and IC contributions to the SED are both parametrized in logiv) — 
log(L) space by a linear curve which transitions to a parabolic curve, para metrized by the 
151 MHz luminosity, Li 5 i MHz ( IFossati et al.lll997l . Il998l ; iDonato et al.ll200ll ). For x = logv, 
xr = log(151 MHz), xx = log(l keV/hp), where h p is Planck's constant, ip(x) = logL(logzz), 
iJ) R = log L 151M hz, and ip x = \ogL lkeV , L[y) can be found by 



where 



and 



i/>(x) = log(10^ s + 10^ c ) 

(1 - a s )(x - x R ) + ijj R x<x trs 
-[(x - x s )/a} 2 + i) sp x > x tTS 



(Al) 
(A2) 



^ _ | (1 - a e )(x - x x ) x<x trc 

' " \ —[(X - X c )/<7] 2 + 4l cp X > X trc 

Fossati et al.l (119971 ) gives the synchrotron slope a s = 0.2, the synchrotron transitional fre- 
quency v t rs = 5 x 10 10 Hz, and the width parameter 



a 



Xtrs X s 

1 - a. 



-,1/2 



(A4) 



while IDonato et al.l (120011 ) gives the inverse Compton peak slope a c = 0.6. 
To ensure continuity of ips at x trs , the parameter ip sp is defined 

-i 2 

*^ s 

i) sp = (1 - a s ){x t rs - Xr) + 



a 



+ if) 



R- 



(A5) 



For ipic to be continuous and differentiable at x trc , ipx must be defined as 

— x 1 2 

- (1 - a c )(x trc - x x ) + ipc P 

a 



ip x = - 



Xfr 



(A6) 



and the transition frequency between the linear and parabolic portion of ipic, v trc , must be 
defined such that 

Xtrc = x c - a 2 (l - a c )/2. (A7) 

In §2.1.31 the synchrotron peak frequency v$ = 10 xs , the ratio between the inverse Compton 
and synchrotron peak frequencies fic/^s — x c — x s , and the ratio between the inverse Comp- 
ton and synchrotron peak luminosities Ljc/L$ = ipcp — ?Psp are set according to observational 
data and individual blazar SED models. 
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Tab le 1. FSR Qs used for dete r minin g L T c / L$. Sources used in t able: F06 



(120061) . T0 7lTavecchio et a 
(bopgh . Gi09 iGiuliani et alH 



to07|), M0 8 Maraschi et a 
2009h. Gh 09lGhisellini et al 



(12009( 1. & W09 IWatanabe et all (l2009h 



Foschini et al. 



(2008), B09 iBianchin et al 



(2009b), V09 IVercellone et al 



Name Source « log(Lic/Ls) 



RBS 315 


T07 


2.0 


S5 0836+71 


T07 


2.0 


J0746.3+2548 


W09 


1.5 


PKS 2149-306 


B09 


1.0 


3C 279 


Gi09 


1.0 


3C 454.3 


V09 


1.0 


2141.2+1730 


M08 


0.0 


0521.7+7918 


M08 


0.0 


1234.9+6651 


M08 


1.0 


1050.9+5418 


M08 


0.5 


0402.0-3613 


M08 


0.0 


0828.7+6601 


M08 


1.0 


1623.4+2712 


M08 


1.0 


1340.7+2859 


M08 


1.0 


0152.4+0424 


M08 


1.5 


0232.5-0414 


M08 


0.0 


SDSS J081009.94+384757 


M08 


1.5 


MG3 J225155+2217 


M08 


3.0 


0048-071 


Gh09 


1.5 


0202-17 


Gh09 


1.5 


0215-015 


Gh09 


0.0 


0528+134 


Gh09 


2.0 


2251-158 


Gh09 


1.0 


0227-369 


Gh09 


2.0 


0454-234 


Gh09 


1.0 


0347-221 


Gh09 


2.0 


0820+560 


Gh09 


1.5 


0917+449 


Gh09 


1.0 


1454-354 


Gh09 


1.0 


1013-054 


Gh09 


1.5 


1502+106 


Gh09 


1.0 
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Table 1 — Continued 



Name 


Source 


« \og{Lic/Ls) 


1329-049 


Gh09 


1.5 


1520-319 


Gh09 


2.0 


1551+130 


Gh09 


1.0 


2052-447 


Gh09 


2.0 


1633+382 


Gh09 


1.5 


2227-088 


Gh09 


1.5 


2023-077 


Gh09 


2.0 


2325+093 


Gh09 


2.0 


PKS 1334-127 


F06 


0.5 



Table 2. HBLs used for de termining Lm/Ls 



Aharonian et al. (2008), C08 ICostamante et al. 



Tagliaferri et all (120081 ) . G09 



Ghiscllini ct al 



and vj cjvs- Sources used in table: A08 
J2005), M0 8lMaraschi et al (l2008h . T08 
f 2009b ). F06 iFoschini et all (hoOfih " MOB 



Massaro et al. 



fl2006h . 



Name 


Source 


« log(Lj C /L s ) 


~ log(i//c/ vs) 


RGB J0152+017 


A08 


0.0 


8.0 


PKS 2155-304 


C08 


0.0 


9.0 


J1456.0+5048 


M08 


0.0 


9.0 


1ES1959+650 


T08 


-1.5 


8.0 


0426-380 


G09 


1.5 


8.0 


0235+164 


F06 


0.0 


8.0 


Mkn 501 


M06 


0.0 


8.0 
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Table 3. Parameters of AGN X-ray Luminosity Functions considered. a in units of h\ 
Mpc -3 . b in units of erg s _1 . c an LADE model. 



Luminosity Function 


AjA% 


logL*/logL b 


7i 


72 


Pi 


P2 


Z*c/Zc 


log^ 


a 


Ueda et al. (2003) 


5.04e-6 


43.94 


0.86 


2.23 


4.23 


-1.5 


1.9 


44.6 


0.335 


La Franca et al. (2005) 


1.21e-6 


44.25 


1.01 


2.38 


4.62 


-1.15 


2.49 


45.74 


0.20 


Silverman et al. (2008) 


6.871e-7 


44.33 


1.10 


2.15 


4.22 


-3.27 


1.89 


44.6 


0.333 


Ebrero et al. (2009) 


4.78e-6 


43.91 


0.96 


2.35 


4.07 


-1.5 


1.9 


44.6 


0.245 


Yencho et al. (2009) 


7.24e-7 


44.40 


0.872 


2.36 


3.61 


-2.83 


2.18 


45.09 


0.208 


Aird et al. (2009) c 


2.95e-5 


44.77 


0.62 


3.01 


6.36 


-0.24 


0.75 




-0.19 



Table 4. f C r needed for luminosity function to match the peak of the X-ray background 
and corresponding CT number density at z = in Mpc -3 for L x > 10 43 erg s -1 



Luminosity Function /ct CT number density at z = (Mpc 3 ) 

Without Blazars With Blazars Without Blazars With Blazars 



Ueda et al. (2003) 


0.5 


0.4 


7.3 x 10~ 6 


4.4 x 10~ 6 


La Franca et al. (2005) 


0.3 


0.2 


2.0 x 10" 6 


1.6 x 10~ 6 


Silverman et al. (2008) 


0.8 


0.8 


1.4 x 10~ 5 


1.2 x 10" 5 


Ebrero et al. (2009) 


0.1 


0.02 


1.1 x HT 6 


3.6 x 10- 7 


Yencho et al. (2009) 


0.8 


0.8 


1.7 x 10~ 5 


1.6 x 10~ 5 


Aird et al. (2009) 


0.7 


0.6 


5.2 x 10~ 6 


4.7 x 10~ 6 
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Fig. 1. — Rest frame radio AGN lumino sity function by Willo tt et al. (2001) model C rela- 
tivistically beamed using the method of lUrry fc Shaferl (119841 ) and lUrry fc Padovanil (119911 ) 
at z — 1 separated into FSRQs (solid line) and BL Lacs (dashed line). 
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log(E) (keV) 

Fig. 2. — Rest frame spectral energy distributions used for FSRQ (solid line) with log zaL„(151 
MHz)=43.0, LBL (dashed line) with logz/L„(151 MHz)=41.5, and HBL (dot-dashed line) 
with logi/L„(151MHz)=40.0. 
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Fig. 3. — AGN and Blazar contribution (solid line) to the X-ray and 7-ray background if 
BL Lac duty cycle is 100%. AGN (dotted line- using Ueda et al. 2003 X-ray luminosity 
function and CT fraction f ct = 0.5), FSRQs (dashed line), and BL Lacs (dot-dashed line). 
The colored data and areas denote measurements from various instruments: blue - ASCA 
GIS (Kushino et al. 2002); magenta - RXTE (Revnivtsev et al. 2003); green - XMM- Newton 
(Lumb et al. 2002); red - BeppoSAX (Vecchi et al. 1999); yellow - ASCA SIS (Gendreau et 
al. 1995); cyan - XMM- Newton (De Luca & Molendi 2004); light yellow - SMM (Watanabe et 
al. 1999); grey data - HEAO-1 (Gruber et al.1999); blue data - INTEGRAL (Churazov et al. 
2007); red data - SWIFT BAT (Ajello et al. 2008); cyan data - COMPTEL (Weidenspointer 
et al. 2000); red triangles - EGRET (Sreekumar et al. 1998); blue triangles - reevaluation 
of EGRET (Strong et al. 2004); green data - renormalization of EGRET based Sreekumar 
et al. 1998 (Stecker et al. 2008). 
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Fig. 4. — AGN and Blazar contribution (solid line) to the X-ray and 7-ray background. 
AGN (dotted line- using Ueda et al. 2003 X-ray luminosity function and CT fraction f ct = 
0.4), FSRQs (dashed line), and BL Lacs (dot-dashed line) with an X-ray duty cycle of 13%. 
Data the same as in Figure [3j 
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Fig. 5. — BL Lac number density counts for 15-55 keV band assuming an X-ray duty cycle 
of 13%. Data shown from Ajello et al. (2009) Figure 12b. 
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Fig. 6. — Number d ensity of CT AGN with Lx > 10 43 erg s -1 as a functio n of redshift for 
the HXLFs given by lUeda et al.l (120031 ) (solid lines) and lEbrero et al.l (120091 ) (dashed lines). 
The black(thin) lines are the CT AGN needed if blazars are not considered. The blue(thick) 
lines are the CT AGN needed if the blazar contribution to the CXB is considered. 



